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We discuss a model for the space-time evolution of ultrarel- 
ativistic heavy-ion collisions which employs relativistic hydro- 
dynamics within one region of the forward light-cone, and mi- 
croscopic transport theory (i.e. UrQMD) in the complement. 
Our initial condition consists of a quark-gluon plasma which 
expands hydrodynamically and hadronizes. After hadroniza- 
tion the solution eventually changes from expansion in local 
equilibrium to free streaming, as determined selfconsistently 
by the interaction rates between the hadrons and the local 
expansion rate. We show that in such a scenario the inverse 
slopes of the m-r-spectra of multiple strange baryons (S, Q) 
are practically unaffected by the purely hadronic stage of the 
reaction, while the flow of p's and A's increases. Moreover, 
we find that the rather "soft" transverse expansion at RHIC 
energies (due to a first-order phase transition) is not washed 
out by strong rescattering in the hadronic stage. The earlier 
kinetic freeze-out as compared to SPS-energies results in sim- 
ilar inverse slopes (of the mr-spectra of the hadrons in the 
final state) at RHIC and SPS energies. 



Ultrarelativistic heavy ion collisions offer the unique 
opportunity to study highly excited QCD-matter in the 
laboratory, and possibly the QCD phase transition to the 
so-called quark-gluon plasma (QGP) at high energy den- 
sity G|. The dynamics of such reactions is commonly 
described either within hydrodynamics or within micro- 
scopic transport models [Q. The fluid-dynamical ap- 
proach is most accurate in very dense systems (where 
mean free paths are small) close to local equilibrium. 
In this model the reaction dynamics is closely linked to 
the equation of state (EoS), which enters directly into 
the equations of motion. It is therefore well suited to 
study systems that undergo phase transitions, e.g. that 
of hadronic matter into a QGP. In particular, one can 
employ an EoS with a first order phase transition, where 
the QGP coexists with hadronic matter within some re- 
gion of energy- and baryon density. 

Hadronic cascades based on binary hadron-hadron col- 
lisions are probably less well suited to study very dense 
systems, and of hadrons coexisting with a QGP. However, 
they are useful for the dilute space-time regions since they 
account for finite relaxation times in the hadron gas |^] 
and the breakup of local kinetic equilibrium. Moreover, 
since each hadron is propagated individually, and its in- 
teractions with other hadrons are described on the basis 



of elementary processes, microscopic transport models of- 
fer the opportunity to calculate the freeze-out conditions 
instead of just putting them in by hand as is usually done 
in fluid-dynamical approaches. In particular, one needs 
not assume that all hadron species decouple on the same 
hypersurface 

The sequencial freeze-out is particularly important to 
understand || the experimental fact that in Pb+Pb col- 
lisions at CERN-SPS energy (^/s = 17 A GeV) the trans- 
verse mass spectra of multiple strange baryons [0J|] are 
"softer" than expected by a linear interpolation of the in- 
verse slopes of protons and deuterons Q: If all hadrons 
flow with the same collective velocity and decouple on the 
same hypersurface, the inverse slopes or average trans- 
verse momenta have to increase with the mass of the 
hadron (at least if one neglects resonance decays). 

We shall discuss here, whether it is possible to under- 
stand these observations by assuming the formation of a 
QGP which expands as an ideal relativistic fluid, going 
through a mixed phase where it coexists with hadronic 
matter. After hadronization is completed (locally in 
space-time), we evolve the system within microscopic 
transport theory (we employ the UrQMD model and 
follow the evolution of the hadrons until they freeze-out. 
We point out that the hadrons are "born" in an expand- 
ing three- volume (i.e. the hadronization hypersurface), 
the expansion rate of that volume being many orders 
of magnitude larger than in the universe at hadroniza- 
tion JlT[ . It will thus depend on the corresponding ele- 
mentary cross-sections and the composition of the system 
which hadrons are able to maintain local kinetic equilib- 
rium, despite the very large expansion rate. We will show 
in particular how the momentum distributions of various 
hadron species are affected by the hadronic stage as com- 
pared to those on the hadronization hypersurface. 

Let us first briefly describe the specific form of hy- 
drodynamics that we shall employ. For a more de- 
tailed discussion see jl2|,[n|. For simplicity, we assume 
boost-invariant longitudinal flow, v z = z/t. This should 
be a reasonable first approximation around midrapidity. 
Cylindrically symmetric transverse expansion is super- 
imposed. For T > Tc — 160 MeV we employ the well- 
known MIT bagmodel equation of state, assuming for 
simplicity an ideal gas of quarks, antiquarks (with masses 
m u = m d — 0, m s — 150 MeV), and gluons. For T <T C 
we employ an ideal hadron gas that includes the complete 
hadronic spectrum up to a mass of 2 GeV. At T = Tc 



we require that both pressures are equal, which fixes the 
bag constant to B — 380 MeV/fm 3 . By construction the 
EoS exhibits a first-order phase transition. 

For collisions at SPS energy we assume that hydrody- 
namic flow sets in on the proper time hyperbola r, = 
1 fm/c. This is a value conventionally assumed in the lit- 
erature, cf. e.g. 0,|l3|]. We further employ a (net) baryon 
rapidity density (at midrapidity) of dNs/dy — 80, as 
obtained by the NA49-collaboration for central Pb+Pb 
reactions Q . The average specific entropy in these colli- 
sions is ~s/~Pb — 45 ± 5 (the bar means averaging over the 
transverse plane) . With this entropy per net baryon most 
measured hadron multiplicity ratios can be described 
within ±20% |lf|. The corresponding initial energy and 
net baryon densities (e, = 6.1 GeV/fm 3 , = 4.5po) are 
assumed to be distributed in the transverse plane accord- 
ing to a so-called "wounded nucleon" distribution with 
transverse radius Rt = 6 fm ]D|] . 

The measured p^-spectra (at midrapidity) of 7r's, p's, 
A's, </>'s, and cfs can be reproduced within this model if 
one assumes an ideal hydrodynamic expansion between 
the t = n hyperbola and the T = 130 MeV isotherm jl3) . 
However, the resulting spectra of S's and fl's are harder 
than those of p's and A's, because the same freeze-out 
hypersurface is employed for all hadrons. 

Due to the higher parton density at midrapidity (as 
compared to collisions at SPS energy), thermalization 
may be reached earlier at RHIC. According to various 
studies M), it may occur within 0.2 — 1 fm. We em- 
ploy Tj = 0.6 fm. The net baryon rapidity density 
and specific entropy at midrapidity in central Au+Au 
at = 200A GeV is predicted by various models of the 
initial evolution to be in the range dNs/dy « 20 — 35, 
s/p B Ri 150-250 0. We will employ dN B /dy = 25 and 
s/p B = 205 (-► €i = 20 GeV/fm 3 , p t = 2.3p ). (The ini- 
tial conditions could of course be fine-tuned once the first 
experimental data are available.) This yields a transverse 
energy on the r = Tj hyperbola of dEx/dy = 1.3 TeV, 
which decreases to 720 GeV on the hadronization hyper- 
surface. We find that dEr/dy is almost conserved during 
the following hadronic evolution, i.e. dissipation compen- 
sates for the expansion work. 

Having specified the initial conditions and the EoS, 
the hydrodynamical solution between the r = t% hyper- 
bola and the hadronization hypersurface is uniquely de- 
termined. For a more detailed discussion of how this 
solution is obtained, cf. Jl3[ . 

The number of hadrons of species i hadronizing at 
space-time rapidity r/, proper time r, and position 
rr (sin (x — 0) j cos (x — <^))j with four-momentum p^ = 
(rax cosh y, px sin x, Pt cos x, m T sinh y) , is given by |Q8[ 
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where = (cosh r), vt sin(x — 4>), vt cos(x — <f>), sinh. 77) 
denotes the fluid four-velocity. £ £ [0, 1] parametrizes the 
hadronization hypersurface^] (counter-clockwise), such 
that Vt(0 and T (C) specify the space-time points on the 
hypersurface. / is either a Bose-Einstein or Fermi-Dirac 
distribution function. We assume that vt, T, and the 
chemical potentials, are independent of 77, such that the 
hadronization hypersurface in the t — z plane (fixed tt) 
is simply a proper-time hyperbola. 

The ensemble of hadrons generated according to 
eq. (0.1) is now taken as input for the microscopic trans- 
port model UrQMD. Thus, after formation each hadron is 
propagated individually along a classical trajectory and 
can scatter stochastically until it eventually does not in- 
teract any more and freezes out. Note that by construc- 
tion the hadron fluid starts from a state of local equi- 
librium, but in a three-space with very rapidly increas- 
ing volume-measure ]Tl] ] . The energy- momentum tensors 
and conserved currents (of UrQMD and hydrodynamics) 
match on the hadronization hypersurface, because the 
local collective flow velocities, as well as the energy and 
net baryon densities, and the pressures, are equal (we 
include the same states in the hadronic part of the fluid- 
dynamical EoS as in UrQMD). 

The UrQMD model is based on the covariant prop- 
agation of all hadrons, excitation of resonances and 
strings and their subsequent decay resp. fragmentation. 
Free cross-sections for hadron-hadron scattering are em- 
ployed. Comparisons of UrQMD calculations to various 
experimental data from SIS to SPS energies, as well as a 
detailed description of the model, are documented else- 
where 0. 

Fig. [jjcompares the my-spectra on the had roni zation 
hypersurface (thin lines), obtained from Eq. (3.1) (plus 
strong resonance decays) , with those at freeze-out (thick 
lines). The (+)-(-) (lj, A, A |o), and 3 @ spectra 
of NA49, and the A, A, S, and ft spectra of WA97 @ 
are also depicted (the latter have been normalized to our 
calculation). Clearly, because of rescattering, the trans- 
verse flow of TV's and A's increases during the hadronic 
stage. Also, slow antinucleons (with m-r < 1.5 GeV) 
are annihilated to a large extent because the annihila- 
tion cross-section increases rapidly at small momenta. 
On the other hand, the spectra of fi's and of S's with 
rriT > 1.6 GeV are practically unaffected by the hadronic 
stage and closely resemble those on the phase boundary. 
This is due to the fact that the scattering rates of 2 and 
f2 in a pion-rich hadron gas are smaller than those of A^'s 
and A's 0. Note that a n and Q can not be coupled to 
any of the known plj B = 1, S = — 3 states with masses 
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1 The switch is actually done slightly behind the hadroniza- 
tion hypersurface, cf. the discussion in fll3|], because the rate 
of expansion of its volume can diverge [pj . 
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< 2.5 GeV. The excitation of a state above this thresh- 
old requires a large (AE ~ 1 GeV) energy transfer from 
the 7r to the SI (in the local rest-frame), which is highly 
suppressed since on the hadronization hypersurface the 
7r-energies are distributed thermally with (E) v ~ 3Tc- 
Finally, elastic collisions give much smaller energy and 
momentum transfers to the S7 than toq , and therefore do 
not change its m^-distribution. 
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FIG. 1. Transverse mass spectra of N, N, A + E°, 
H° + H — , and Q~ at y c .m. — in central Pb+Pb collisions at 
y/a = 17.4 GeV (lower panel) and Au+Au at Vs = 200A GeV 
(upper panel) . Thin lines: on the hadronization hypersurface; 
thick lines: at freeze-out; diamonds: NA49 data; squares: 
WA97 data. 
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FIG. 2. Distribution of the number of interactions that the 
final-state particles suffer after being hadronized. 

Fig. |^ displays the inverse slopes T* obtained by 
a fit of dNij 'cPrnxdy to exp(— mr/T*) in the range 
rriT — mi < 1 GeV (for clarity, the calculated points have 
been slightly displaced to the left or right of the true 
hadron mass). The trend of the SPS data, namely the 
"softer" spectra of 3's and Si's as compared to a linear 
T*(m) relation, is reproduced reasonably well. This is 
in contrast to "pure" hydrodynamics with kinetic freeze- 
out on a common hypersurface (e.g. the T = 130 MeV 
isotherm), where the stiffness of the spectra increases 
with mass, cf. Fig. || and also refs. || |l!j , |3| . Resonance 
decays are not included in the hydrodynamic spectra on 
the T — 130 MeV isotherm. 
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As shown in Fig. g, on average the baryons which fi- 
nally emerge as S's and Si's suffer even less interactions 
than the final-state p's and A's. Thus, within the model 
presented here, these particles are emitted directly from 
the hadronization hypersurface without further rescatter- 
ing in the hadronic stage. The hadron gas emerging from 
the hadronization of the QGP (in these high-energy re- 
actions) is almost "transparent" for the multiple strange 
baryons, especially because of the large expansion rate. 
In our calculation, the flow seen in the Sl-spectra at SPS 
and RHIC is fully accounted for by the expansion pro- 
ceeding hadronization (for entropy and baryon density 
as discussed above). The space-time domains of freeze- 
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FIG. 3. Inverse slopes of the m-r-spectra of tv, K, p, A+E°, 
H° + H~, and Q~ at y c . m . = 0, tut — rm < 1 GeV. 

When going from SPS to RHIC energy, the model dis- 
cussed here generally yields only a slight increase of the 
inverse slopes, although the specific entropy is larger by 
a factor of 4-5 ! The reason for this behavior is the first- 
order phase transition that leads to smaller transverse 
flow velocities than an ideal gas |Q (but to a larger ex- 
pansion rate of the hadronization volume jllj !). For 



our initial conditions, the collective transverse flow ve- 
locity on the hadronization hypersurface increases only 
from « 0.3 (Pb+Pb at SPS) to w 0.35 (Au+Au at 
RHIC) As can be seen from the present calculation, 
this is not counterbalanced by increased rescattering in 
the purely hadronic stage (compare to the inverse slopes 
obtained from "pure" hydrodynamics with freeze-out on 
the T = 130 MeV isotherm !). 

In summary, we have employed relativistic hydrody- 
namics to describe the evolution and hadronization of a 
hypothetical quark-gluon fluid at CERN-SPS and BNL- 
RHIC energies. The produced hadrons are propagated 
within a microscopic transport model (UrQMD). Inter- 
actions within the hadron gas increase the collective flow 
beyond that present at hadronization, and reduce the 
temperature below the QCD phase transition temper- 
ature (we assume Tq — 160 MeV). As an exception, 
we find that multiple strange baryons practically do not 
rescatter within the hadron gas. Their mT-spectra are 
thus determined by the conditions on the hadronization 
hypersurface, i.e. Tc and the collective flow created prior 
to hadronization. Their spectra therefore are less sen- 
sitive to the confined phase, T < Tq, but are closely 
related to the EoS of the QGP and the phase transition 
temperature Tc- At RHIC energy, for all hadrons T* 
or (pt) are considerably smaller than predicted by pure 
hydrodynamics with freeze-out on the T = 130 MeV 
isotherm. Thus, kinetic freeze-out occurs closer to the 
phase boundary than at SPS energy; this is due to the 
different chemical composition of the central region (more 
mesons, less baryons) and due to the larger expansion 
rate of the hadronization volume. 
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